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Abstract 

The periodic motion of a classical point particle in a one-dimensional double¬ 
well potential acquires a surprising degree of complexity if friction is added. Finite 
uncertainty in the initial state can make it impossible to predict in which of the two 
wells the particle will finally settle. For two models of friction, we exhibit the structure 
of the basins of attraction in phase space which causes the final-state sensitivity. 
Adding friction to an integrable system with more than one stable equilibrium emerges 
as a possible “route to chaos” whenever initial conditions can be specified with finite 
accuracy only. 


1 Introduction 

The difficulty to reliably predict the behaviour of a classical systems is usually related 
to the existence of fractal structures in the mathematical model describing the system. 
Conservative non-integrable systems such as three interacting planetary bodies [ij and 
chaotic dissipative systems such as Lorenz’s model of the atmosphere [2] provide two 
well-known cases in point. Consider adding a third body to the integrable system of two 
planetary bodies which interact through gravitation. The KAM theorem [$] describes how 
the original foliation of the system’s phase-space into tori is being replaced gradually by a 
highly intricate mixed phase space. Finite balls of initial conditions will no longer contain 
trajectories on tori only but also others which separate at an exponential rate. The non¬ 
linearity present in Lorenz’s model gives rise to a strange attractor [JJ. Its properties 
dominate the long-term evolution of the system since trajectories with neighboring initial 
conditions are likely to visit rather different regions of phase space at comparable later 
times. This phenomenon has been called final state sensitivity [5]. 

An actual macroscopic physical system, however, cannot exhibit fractal structures in 
a strict sense: the system would need to match its mathematical description on arbitrary 
length scales |6j but classical models break down on the molecular or atomic level. Exper¬ 
imentally observed structures may be highly intricate over many - but not all - orders of 
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magnitude. Nevertheless, finite intricacies are sufficient to make reliable long-term predic¬ 
tions impossible when combined with initial states which are known only approximately. 

Adding friction to conservative non-integrable model systems washes out fractal struc¬ 
tures on the finest scales but remnants may continue to exist. This has been shown for a 
spherical pendulum with three stable equilibrium positions, in the presence of gravity |7j. 
When friction is added, basins of attraction with highly intricate boundaries emerge. If 
the initial state is not known exactly the fine structure of the basin boundaries - in spite 
of not being fractal - already prevents any reliable prediction of the equilibrium position 
near which the pendulum will come to rest. 

The purpose of this paper is to identify a different mechanism which also results in an 
unpredictable final state if initial conditions are known only with finite accuracy. Starting 
from an integrable system with multiple stable equilibria, we will show that the addition 
of friction may create basins of attraction with intricate boundaries, leading to a situation 
which resembles the one of the pendulum just described. The main difference is that this 
mechanism does not rely on pre-existing fractal structures: instead, incorporating friction 
causes intricate structures to emerge in an originally integrable system. 

Sec. [2]of this paper provides an initial, qualitative explanation of why adding friction to 
an integrable system with two or more stable equilibria may cause final-state sensitivity. 
Then, in Sec. [3j we investigate the structure of basins of attraction generated by two 
different types of friction when acting on a particle moving in a piece-wise constant double¬ 
well potential. Finally, we summarize and discuss our results in Sec. [IJ 

2 Final-state sensitivity in a double-well potential 

To illustrate how friction creates final state sensitivity, we consider a classical particle 
moving along a straight line in the presence of a symmetric double-well potential. The 
system is described by the Hamiltonian function 

H(p,q) = ^ + W(q), p,q eR, (1) 

where q and p denote position and the momentum of the particle, respectively. The minima 
of the potential W ( q ) are located at q = q± , separated by a barrier of height Wo = W (0) 
which defines the critical energy, E c = Wo. The phase-space diagram of the system is 
shown in Fig. [l] displaying the familiar types of trajectories. The minima L and R of the 
potential are stable fixed points each surrounded by periodic orbits with energies E not 
exceeding the critical value, 0 < E < E c = Wo- For E = E c , the particle may rest at 
the unstable fixed point at q = 0, or travel on one of the two separatrices connected to it. 
The trajectories with energy above the critical value, E > E c are periodic, encircling both 
minima on a single round trip. With a single degree of freedom and the energy H{p , q) as 
a conserved quantity, the system is integrable, leading to the global foliation of its phase 
space into one-dimensional tori. 

Adding friction will modify all trajectories except when then particle initially rests at 
one of the three fixed points. If located on a separatrix or on any periodic trajectory with 
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Figure 1: Phase-space structure of a double-well potential W(q ) with minima L and R, 
located within wells £ and 1Z (shaded areas), respectively; a particle with neighbouring 
initial conditions z, z! (full and dashed lines, respectively) may end up in different wells 
when subjected to periodic, dissipative “kicks” which reduce its momentum and thus its 
energy. 

energy less than E c , dissipation will cause the particle to “spiral” into either the left or 
the right minimum of the potential W(q), depending on its original position relative to 
the origin, q = 0. The particle cannot escape from a well once it has been trapped, and 
the fixed points L and R turn into attractors. 

The destiny of a particle with initial energy E > E c , however, it is not immediately 
obvious since it may end up in either well. Friction will inevitably “draw” the particle 
towards the location of the separatrices of the unperturbed system. At some time, the 
energy of the particle will drop below the critical value E c . The position of the particle 
relative to the origin at the time of the drop will determine whether it becomes trapped 
in the left or in the right well. 

For simplicity, let us assume that friction acts at discrete times only, repeatedly reduc¬ 
ing the momentum of the particle by a constant factor. Suppose that for initial conditions 
z = (Pi q ) T , the particle will - after a possibly long time - settle in the right well as 
illustrated in Fig. [lj Intuitively, a slightly smaller initial momentum (see zl in the figure) 
could cause the particle to negotiate the barrier one less time and to settle in the left well 
instead. The slight change in the initial condition has thus altered the long-term behaviour 
of the system. Therefore, the finally state of a particle may become unpredictable from a 
practical point of view, i.e. whenever its initial conditions are known to lie within a small 
but finite volume of phase space only. 

More formally, the non-Hamiltonian equations of motion map an initial state z(to) to 
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a new value z (t) at time t, 


z(to) ^ z(i), 


z = 


P 

q 


( 2 ) 


leading to a decrease of the energy defined in |T]): Eq —>• E < Eq. To ascertain whether 
a particle with initial state z(to) ends up near L or near R, one needs to determine the 
earliest time t such that its energy E falls below the critical value, 

E <E C . (3) 

Repeating this calculation for all initial conditions will divide the phase space into two 
disjoint sets known as basins of attraction which encode whether the particle ends up in 
well L or 1Z. Let us investigate the structure of their boundaries for two models of friction, 
using a particularly simple double-well potential. 


3 Piece-wise constant double-well potential with friction 

The double-well potential considered here is based on a “particle in a box” defined by two 
infinitely high potential walls at q = ±1 which restrict motion to a line segment of length 
2£. The particle bounces off the walls elastically resulting in an instantaneous reversal of 
its momentum: p —> —p] its position q = ±1 remains unchanged when hitting a wall. A 
piece-wise constant potential, 


W(q) = 



\Q\<e, 

£ < \q\ <1, 


(4) 


models the smooth double-well. For simplicity, we take an arbitrarily thin potential bar¬ 
rier, corresponding to e —> 0. The only impact of this “infinitesimal” barrier is to confine 
the particle in a well once its energy drops below the critical value Wo, thus creating the 
wells C and 7Z. Two continuous sets of potential minima exist because the bottom of the 
potential is flat. 

A widespread method to investigate non-integrable systems is to start with an inte¬ 
grate system and add a perturbation, be it a time-independent potential term as in the 
KAM theorem or a deterministic time-dependent force |8]. To support our claim that fric¬ 
tion generically causes final state sensitivity, we will model it in two different ways which 
are inspired by these approaches. In the first case, the elastic collisions of the particle 


with the boundary walls are made inelastic (cf. Sec. 3.1) while an impulsive friction force 
is applied periodically in the second case (cf. Sec. 3.2). The first model depends on a 
single parameter only, the coefficient of restitution. The second model depends on two 
parameters, the frequency and the strength of the dissipative “kick.” 


3.1 Inelastic collisions 

The motion of the particle in the piece-wise constant double well © consists of free 
motion between the walls interspersed with momentum-reversing elastic collision at the 
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walls. The dynamics changes fundamentally upon replacing the elastic collisions at the 
walls by inelastic ones, characterized by a coefficient of restitution, r E (0,1): 


z i —> z = R z , R = 



for q = ±£. 


(5) 


This minor change turns the conservative system into a dissipative one and - as we will 
see - is sufficient to create an embryonic form of final-state sensitivity. 

The long-term dynamics of the particle does not depend on its initial position: all 
particles with fixed momentum p$ but arbitrary position qo E (—£,£) will experience the 
same amount of friction, only to end up in same well. Thus, let us assume that the particle 
starts out with positive initial momentum po > p c , beings located at qo = £-, i.e. just to 
left of the right wall. Then, the initial state zq = (p q,£~) at time to evolves according to 

z (i„) = R n • z 0 = f P ° V n E No , (6) 


with the times t n being defined by particle returning to its initial position £_. Monitoring 
the value of its momentum at the walls is sufficient to determine the well which will trap 
the particle. The particle will be trapped in well 7Z, for example, if its last collision at the 
right wall makes its energy drop below the critical value E = E c due to p H > ( —r)p. 

For positive initial momentum po the particle will hit the right wall first. The well to 
finally trap the particle is determined by the number of collisions ti{Eq) before it drops 
below E c . Denoting the energy of the particle after n collisions by E n , we need to find the 
number n(E$) such that the energy of the particle falls below E c for the first time, 


n(Eo) = min {E n < E c } , with E$> E c . 

n£ N 


Using Eq. ([h]) the number n(Eo) is easily found to be 


n(E 0 ) 


1 In (E c /Eq) 


In (pc/po) 

2 lnr 


lnr 


(7) 

( 8 ) 


where the initial moment defines the initial energy, Eq = Pq/2?7t,, and \x\ is the ceiling 
function extracting the smallest integer greater or equal to the number x. If n is odd 
(even), a particle with positive momentum po will end up in the well on the right (left). 
The basins of attraction for the wells L and 1Z are given by alternating horizontal bands 
in phase space shown in Fig. [2j The widths of the bands decrease with decreasing friction 
(and they increase with energy E which the figure does not show due to the limited 
momentum range). If the initial conditions (po,qo) of a particle are known exactly , then 
the deterministic dynamics leads to a unique and well-defined final state which can be 
predicted with certainty. However, limited precision of the initial conditions may results 
in a genuine indeterminacy of the final state. Assume that the initial state of the particle 
is only known to lie inside a rectangle with sides A q > 0 and A p > 0, centered about 
the point zq. Trajectories with initial momenta po and p ' 0 = rpo are bound to end up in 
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Figure 2: Basins of attraction for the phase-space region (1 < p/p c < 1.2, 0 < q/i < 0.1) 
of a particle of unit mass in a double-well, with friction arising from inelastic collisions at 
the boundary walls: initial conditions located in dark (light) regions will end up in the 
right (left) well of the potential. The vertical bars correspond to different values of the 
coefficient of restitution: r = 0.9, r = 0.99, r = 0.999 (left to right). Smaller values of 
friction lead to increasing “complexity” of the basin boundaries, in the sense of producing 
narrower bands. 

different wells. Thus, if the inaccuracy in momentum exceeds this value, 

Ap> p 0 - p 0 , (9) 

the uncertainty rectangle will cut across at least two adjacent basins of attraction. In 
other words, given the initial momentum po and any finite uncertainty A p about it, the 
prediction of the final state becomes impossible for a coefficient of restitution in the interval 

A p , . 

1- -<r< 1, (10) 

Po 

since the rectangle with sides A q and A p will contain trajectories destined for the wells C 
and 7 Z. We conclude that sufficiently weak inelasticity prevents the reliable prediction of 
the final state. In this well-defined sense, adding friction to an integrable system provides 
a mechanism which prevents accurate long-term predictions. 

3.2 Periodic damping 

Now we turn to a model where friction is caused by a periodic, dissipative force which 
acts during a short time interval only. It will be convenient to consider the limit of an 
instantaneous action which multiplies the momentum of the particle by a constant factor 
7 £ (0,1) at times = kT, with k £ N, and a free parameter T. This approach is 
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analogous to periodically kicking a system with a deterministic force which, for a particle 
moving freely on a ring known as a “rotor,” produces deterministically chaotic motion 
(3]. Since our model depends on two independent parameters, 7 and T, we expect more 
complicated basins of attraction compared to the model with inelastic reflections. 

To construct the basins of attraction of the wells C and 1Z , we need to determine when, 
for arbitrary initial conditions {po,Qo) T , the energy of the particle falls below the critical 
value for the first time We then record whether, at that moment of time, it is located to 
the left or to the right of the origin, i.e. within C or 1Z. For simplicity, the particle is 
assumed to begin its journey at time t = 0 + , i.e. just after t = 0 , with positive momentum 
Po > p c and arbitrary initial position go £ (—£,£)■ 

The particle moves freely during intervals of length T, with perfectly elastic collisions 
occurring at the boundary walls which only change the sign of its momentum. An ex¬ 
pression for its time evolution in closed form can be found if we “unfold” the trajectory 
by imagining identical copies of the double-well to be arranged along the position axis. 
Instead of being reflected at the right wall located at q = £, the particle enters the next 
double well, which occupies the range (£,3£), and continues to move to the right, etc. In 
this setting, the momentum does not change its sign when the particle moves from one 
double well to the adjacent one. The sign of its momentum in the original double well 
is negative (or positive) if the particle has reached the s th copy of the double well, with 
s E N being odd (or even). 

To determine the dynamics of the system over one period of length T, we combine the 
free motion with the dissipative kicks: 

1 . during the motion of the particle from t = 0 + to just before the first kick at time 
t = T, its phase-space coordinates are given by 


z (t) 


Po 

.qo+Pot/m, 


F(t)-z 0 , 


F (i) 




i€(0+T-), (11) 


where q E (0, 00 ) due to the unfolding; 

2. the dissipative kick at time T reduces the momentum of the particle by the factor 

7 6(0,1), 


z(T+) 


fr/p{T )\ 

Ucny 


D-z(T-), 



te(T-,T + ). (12) 


To obtain the actual position and momentum of the particle inside the box at time t, we 
map (or “fold back”) the expression F(t) • z to the interval q E (—£,£), by writing 


z(t) 


H s(t) Po ' 

[(go + Pot/m ) mod 2£] - £, 


te (0+,T+), 


(13) 


where the value of the integer s(t) is determined by writing q + pt/m = q(t) + 2 £s, with 
q(t) E (—£,+£). The momentum p changes sign whenever the “unfolded” coordinate passes 
through the values £, 3£, 5£,... 
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The time evolution of the initial state zo from time t = 0 + to t = Tjf = ( kT ) + , i.e. 
just after the kick with label k, follows from concatenating Eqs. © and (| 12 [) k times, 


! ®=< D ' F ( r »‘^U", 


• Z 0 = 


( 


7 * 0 

\vk{l)T/m 1, 


■ zo , 


where 


°k{l) = 


1 -7 fc 
1-7 


I’GfJ. 


(14) 


(15) 


In analogy to Eq. (13), the “true” coordinates of the particle inside the box are obtained 
as 

. ( -r ( : HV ... „v d«) 


J(<?o + Vkh)PoT/m) mod 2 i\-i) 

assuming that, after k kicks, the energy E k = p\/2m of the particle has not yet dropped 
below the critical value E c . 

We are now in the position to determine which initial conditions zq will send the 


particle to the left and the right well, respectively. Using Eq. (16), we first determine 


the smallest value of k which reduces the energy of the particle below the critical value, 
E k < E c , or 


k r = 


'1 \n(E c /E o y 


'ln(p c /poY 

2 In 7 


In 7 


k c £ N, 


(17) 


assuming, of course, that po > p c - This relation structurally resembles the result (| 8 j) , with 
the number k c of dissipative kicks playing the role of the number of inelastic collisions 


n c . The sign of the position coordinate after k c kicks, q(T, j~), follows from Eq. (16) and 


determines whether the particle is trapped in C or 77. The explicit dependence of z (Tjf) 
on the initial position implies that changes in go may also produce different final states, 
in contrast to the model studied in Sec. E3J 

Fig. i which has been generated numerically on the basis of Eq. ( |16| ), illustrates 
these conclusions. The first vertical bar visualizes the basins of attraction associated with 
the wells C and 77, respectively. The expected dependence on both initial momentum and 
position becomes clearly visible in the magnifications which also reveal that the boundaries 
of the apparently irregular basins of attractions are not fractal. 


The boundaries of the basins can be found directly from Eq. (16): all initial conditions 


( po, go) mapped to a fixed value of position at time kT are located on lines of the form 

in 1 m 

Po(qo) = - 7—= g 0 + const ~--—g 0 + const, 7 < 1, (18) 


using a kip/) — ( 1 /fc) + 0(1 — 7 ), which holds for weak damping, i.e. for 7 approaching the 
value one from below. Consequently, the boundaries of the basins of attraction are straight 
lines in phase space just as for the model with inelastic reflections off the wall. The lines 
are no longer horizontal but their slope approaches the value zero if a large number of 
kicks is required for the particle to settle in a well. 

Assume once again that the initial conditions of the particle can be prepared with finite 
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Figure 3: Basins of attraction for the phase-space region (1 < p/p c < 11, 0 < q/i < 0.2) 
of a particle with unit mass in a double-well, with friction arising from periodic dissipative 
kicks at times kT, k E N, with t = 1, T = 100s and 7 = 0.99: initial conditions located in 
dark (light) regions will end up in the right (left) well of the potential. Each of the three 
vertical bars on the right magnifies a horizontal strip of the bar to its left by a factor of 
ten. The basin boundaries clearly exhibit both a momentum and a position dependence. 
Each bar results from iterating 501 x 501 regularly spaced initial conditions inside the area 
shown. 


precision only, i.e. they lie inside a phase-space rectangle with area AqAp > 0 and center 
zq. For any finite imprecision one can always find a damping strength 7 such that at least 
one basin boundary crosses the rectangle; this is sufficient to prevent the prediction of the 
well to finally trap the particle. For large initial momenta jjq, the reasoning behind the 


derivation of the inequality ( 10 ) also applies here since the strips constituting the basins 


of attraction will, typically, have almost horizontal boundaries. Thus, for any initial 
conditions (po,qo) and finite uncertainties, damping strengths within the interval 


A p 

1 -- < 7 < 1 

Po 


(19) 


correspond to a situation with an unpredictable final state. Occasionally, the uncertainty 
rectangle with sides A q and A p may cover an area where a slight change in initial position 
causes the particle to reach different wells, which only increases the final state sensitivity. 


4 Summary and conclusions 

We have shown that adding friction to an integrable one-dimensional double-potential well 
causes its dynamics to exhibit a rudimentary form of final-state sensitivity. For simplicity, 
the well has been modeled as a “box” divided into two regions by a thin wall. A particle 
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has been subjected to two types of dissipative forces which, by reducing its initial energy, 
cause the particle to necessarily settle in one the wells after a finite, possibly long time. 
The main result of our study is that adding friction to an integrable system produces 
basins of attraction with finely structured boundaries. 

If the particle collides inelastically with the confining walls the resulting basins foliate 
the phase space of the system into horizontal layers of variable width which get narrower 
for decreasing friction. Any ball of initial conditions which extends beyond more than one 
band prevents us from predicting with certainty the well in which the particle will finally 
settle. Periodic dissipative kicks create basins of attraction with slightly more intricate 
boundaries, due to their additional position dependence. Since the particle must settle in 
a well after finite time the observed structures cannot be fractal. In practice, however, it 
is crucial whether the initial conditions can be specified with sufficient accuracy to avoid 
a spread across basins which send the particle to different final states. 

These model systems demonstrate that adding friction to an integrable system with 
multiple stable equilibria can have a fundamental impact on long-term predictability. The 
motion is not “deterministically random” which would require fractal phase-space struc¬ 
tures. However, if the accuracy of the initial conditions falls below a specific threshold, 
the final state of the system cannot be predicted reliably. Experimentally, the precision 
required for a reliable long-term prediction may well be out of reach. 

We expect our conclusions to be structurally stable in the sense that they should not 
depend on the model of friction used. Any dissipative mechanism will, firstly, contract 
all initial conditions into a small phase-space region which is energetically just above the 
barrier of the double well; secondly, the energy of the particle will drop below E c in a 
way which depends sensitively on the initial conditions. Continuous Stokes friction, for 
example, is thus likely to generate similar basins of attraction. 

To systematically study the creation of basins of attraction with intricate boundaries 
in more general, smooth potential wells, we suggest to exploit the existence of action-angle 
variables (/, p>) in integrable systems. The energy E represents a convenient starting point 
to study the impact of friction forces, when expressed as a function of the action /, 

1(E) = j p(q,E)dq. (20) 

Suitable perturbations are easily added to the new form of the Hamiltonian, H = H(I ), 
once position q and momentum p have been mapped to (I, <p) by means of a canonical 
transformation. 

Finally, we highlight a natural application of the effective unpredictability of a final 
state due to friction, given sufficiently imprecise initial conditions. It arises upon intro¬ 
ducing a larger number of identical potential wells arranged on a ring, (37 or 38, say), 
mimicking a one-dimensional roulette wheel. Including periodic dissipative kicks provides 
a surprisingly simple explanation why a finite spread in initial momenta and positions is 
sufficient to generate random outcomes, the working hypothesis underlying any gambling. 
This approach should be contrasted with models of an actual roulette wheel where unpre¬ 
dictable trajectories arise through a multitude of effects: motion in a bent annulus-shaped 
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region embedded into three dimensions, the presence of gravity, rolling resistance and a 

(presumably) non-integrable time-dependent potential. 
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